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ABSTRACT 

We conducted a deep narrowband NB973 (FWHM = 200A centered at 9755A) survey of z = 7 
Lya emitters (LAEs) in the Subaru/ XMM-Newton Deep Survey Field, using the fully depleted CCDs 
newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z = 7 Lya at 
~ 1/mi than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey 
in the Subaru Deep Field that led to the discovery of a z = 6.96 LAE, we detected three probable 
z = 7 LAE candidates. Even if all the candidates arc real, the Lya luminosity function (LF) at 
z = 7 shows a significant deficit from the LF at z = 5.7 determined by previous surveys. The LAE 
number and Lya luminosity densities at z — 7 is ~ 7.7-54% and ~ 5.5-39% of those at z = 5.7 to 
the Lya line luminosity limit of L(Lya) > 9.2 x 10 42 erg s _1 . This could be due to evolution of the 
LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly 
evolves from z = 5.7 to 7. However, even after correcting for this effect of galaxy evolution on the 
decrease in LAE number density, the z — 7 Lya LF still shows a deficit from z = 5.7 LF. This might 
reflect the attenuation of Lya emission by neutral hydrogen remaining at the epoch of reionization 
and suggests that reionization of the universe might not be complete yet at z = 7. If we attribute 
the density deficit to reionization, the intergalactic medium (IGM) transmission for Lya photons at 
z = 7 would be 0.4 < T^^ 1 < 1, supporting the possible higher neutral fraction at the earlier epochs 
at z > 6 suggested by the previous surveys of z = 5.7-7 LAEs, z ^ 6 quasars and z > 6 gamma-ray 
bursts. 

Subject headings: cosmology: observations — early universe — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 



The recent progress on the observations of high red- 
shift galaxies are really remarkable and has been reveal- 
ing how galaxies evolved and contributed to reionization 
Optical narrowband surveys hav e discovered Lya emit 
ters (LAEs) up to z — 6-7 (e.g.. |Shimasaku et al.l l200 



Taniguchi et all 120051: iKashikawa et all 120061: live et al.l 
200a IQuchi et all 120081: IQta et all 120081 ). A few pio- 



neering attempts to search for z ~ 9 LAEs with near- 
infrared narrowband were made and ended with null 
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detections (IWillis fc Courbinl l200l ICubv et al.1 120071) . 
Hight redshift LAEs can be a probe of reionization 
since observed Lya line luminosity function (Lya LF) 
is expected to decline beyond z ~ 6, where reioniza- 
tion is thought to have completed bas ed on the quasar 
and 7-ray burst (GRB) observations (|Fan et al.l 120061 : 
iTotani et al.l 120061 : iGreiner et al.1 120091 ) , as the increas- 
ing fraction of IGM neutral hydrog en absorbs or scat- 
ters t he Lya photons from LAEs (lHaiman fc Spaans I 
[19991 IRhoads fc Malhotra I [200lt IHu et all 120021) . Ear- 
lier works showed that Lya LF does not change much 
beyond z ~ 6 by comapring LFs at z = 5.7 and 6.5 
and concl uded that the un i yerse was largely ionized a t 
z = 6.5 (jStern et al.l [20051: [Malhotra fc Rhoadsl [200l . 
Also, some cosmolog ical simulations of reio nization sup- 
port this result (e.g.. lHaiman fc Cenll2005l) . These stud- 
ies used the samples compiled from various literatures 
with different selection criteria, survey areas and depths 
and suggested that the larger deeper uniform samples 
are necessary to improve the accura cy and statistics 
and d raw t he firmer conclusi o n. La ter, [Shimasak u et all 
( 2006) and Kashi kawa et al.l (|2006l ) constructed a large 
homogeneous samples of z — 5.7 and 6.6 selected from 
the consistent criteria, the same and large survey area 
and comparable depth to find the significant decline of 
Lya LF from z — 5.7 to 6.6 and suggested possible in- 
crease in neutral fraction at z > 6 and the end of reion- 
ization at z ~ 6. For even earlier epoch, iHibon et all 
(2010) recently performed a near-infrared narrowband 
survey and detected seven z — 7.7 LAE candidates. They 
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concluded that Lya LF does not evolve from z — 6.6 to 
7.7 if two brightest candidates in their sample are real 
but otherwise it does evolve, and emphasized that sam- 
ple contamination, if any, easily affects the results and 
thus spectroscopic comfirmation is important. 

Meanwhile, broadband dropout techniques have de- 
tected Lyman break galaxies (LBGs) up to z ~ 6 



10 fe. g..lStanwav et al.M2007t IBouwens et al.ll 2005. 
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2001120101: IHenry et all 120071 120081: IBunker et al.lL 
Hickey et al.l 120091 : lOesch et all 12003 . 120101: lOuchi et al.l 
2009bJ). Their UV continuum LF (UV LF) implies rapid 



buildup of luminous galaxies from z > 7 to z < 6 (e.g., 
iBouwens fc Illingwo rth 2006). However, the faint-end 
slope of z > 6 LBG UV LF and the escape fraction of 
ionizing photons are still highly uncertain, and galaxy 
contribution to reionization has not been tightly con- 
strained yet (IBunker et al.ll2004l:IYan fc Windhorstll2004; 
IBouwens et al.ll2006ll2008t IHenry et al.H2009H . ~ 

Recently, gravitationally lensing galaxy clusters have 
been used to search for magnified extremely faint galax- 
ies at z > 7 and to probe faint end slopes of Lya 
and UV LFs. Three types of ambitious cha llenges have 
been m ade lately. Using the Keck Telescope, IStark et al.l 
(2007) conducted near-infrared blind spectroscopic sur- 
vey of z ~ 8.5-10.4 LAEs by placing long slits to the 
critical lines, locations of the highest magnification, of 
9 lensing clusters, and found six candidates. On the 
other hand, z > 7 LBGs were surveyed over several lens- 
ing clusters using deep optical-infrared imagings of Hub- 
ble Space Telescope (HST), Spitzer and VLT, covering 
most of the magnification regions including critical lines 
(IBouwens et alj|2009t [Bradley et alll2008t IRichard et al.l 



120081 ). IRichard let all 1)20081 ) detected 10 z ~ 7-8 z- 
dropo uts a nd 2z~ 9-10 J -drop outs. Both IStark et al.l 
(2007!) and IRichard et all (120081 ) found that the num- 
ber densities of very faint LAEs and LBGs at z > 7 
are quite high, and they can supply the necessary pho- 
tons to sustain reionization. However, their faintness has 
not yet enabled successful spectroscopic ide ntification of 
them, and they still stay candidates. Finally. fWillis et"aLl 
(2008) conducted a narrowband survey of z ~ 9 LAEs 
targeting 3 lensing clusters and ended with no detection 
of candidates. 

Most of the insights to the galaxy evolution and reion- 
ization at z > 7 are based on candidate galaxies. We 
do not know how many of them are real, and this frac- 
tion easily affects the analysis results. Hence, it is im- 
portant to construct samples of low contamination and 
of spectroscopically confirmed galaxies at z > 7. To 
this end, we conducted a narrowband NB973 (A c = 
9755A and AA = 200l corresponding to 6.94 < z < 7.11 
Lya emission) imaging of the Subsru Deep Field (SDF, 
876arcmin 2 ; [Ka shikaw a~et al.l 1200*41 ) with Suprime-Cam 
(jMivazaki et al.ll2002D on the Subaru Telescope and suc- 
cessfully identifi ed a z = 6.96 LAE by spectroscopy 
(jive et al.|[2006D . The z ~ 7 LAE number density was 
only 17% of that at z = 6.6 with an estimate of lower in- 
tergalactic medium (IGM) transmission for Lya photons 
at higher redshifts; T[f™ ~ 0.62-0.78 at z = 6.6 and 
T Ly<f ~ 0.40-0.64 at z = 7 (iOta et al.M20"08h . This series 
of density declines and decrease in at z > 6 could 

be another evidence that reionization ended at z ~ 6. 
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Fig. 1. — Filter transmission of the Suprime-Cam broadbands 
BVRi'z' and narrowbands NB816, NB921 and NB973 (thin solid 
curves) used for our photometry. The OH night skylines are also 
overplotted with thin dotted curve. The short-dashed curve is the 
atmospheric transmission. The long-dashed and thick solid curves 
show the quantum efficiency of the old CCDs and the new fully 
depleted CCDs. The new CCDs are about twice more sensitive to 
the 2 = 7 Lya at ~ 9730A. 



However, our conclusion from the z = 7 LAE was naive 
because of three weaknesses. First of all, the survey limit 
was as shallow as the Lya luminosity of L(Lya) > 10 43 
erg s _1 , mostly due to the low sensitivity of detector at 
the wavelength ~ 1/im. Hence, our study was limited 
to the very bright end of LF. Secondly, we could detect 
only one z = 7 LAE and thus the sample was statisti- 
cally very small. Of course, this might be really due to 
the significant amount of absorption of Lya photons by 
neutral hydrogen that might remain at z = 7, but we can- 
not completely deny the possibility of the small number 
detection due to the shallow survey depth and the low 
detector sensitivity. Finally, we observed only one sky 
field, the SDF. If we observe any different sky fields, we 
might obtain different results due to the possible cosmic 
variance of galaxy evolution and reionization state. To 
discuss galaxy evolution and reionization more robustly, 
we need a larger z — 7 LAE sample from a deeper survey 
than the SDF over another comparably large sky area. 

The chance to overcome these weaknesses had come 
in 2008. The CCDs of Suprime-Cam were replaced 
with new CCDs that adopt a much larger depletion 
layer, allowing them to detect more photons at redder 
wavelength. The CCDs have twice better sensitivity 
to the z = 7 Lya emission at ~ 1/im than the pre- 
vious CCDs. Inspired with this, we have performed a 
new deeper NB973 filter survey of z = 7 LAEs, tar- 
geting a different but comparably large area sky field, 
the Suba.ru/XMM-N ewton Deep Survey Field (SXDS; 
iFurusawa et al.l 120081 ) . We have detected more candi- 
date z — 7 LAEs to 0.5 mag deeper limit than our previ- 
ous survey in the SDF. In the current paper, we explore 
the photometric properties of them and obtain any im- 
plication for reionization based on the deeper Lya LF 
considering the cosmic variance. The paper is organized 
as follows. In §2, we describe the properties of the new 
CCDs and imaging observation of the SXDS. Then, the 
z = 7 LAE candidates are selected and their photomet- 
ric properties are analyzed in §3. In §4, we derive the 
deeper Lya LF based on the newly detected LAE can- 
didates, compare it with LFs from other recent studies 
and discuss its implication for galaxy evolution and cos- 
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mic reionization. We conclude and summarize our results 
in §5. Throughout, we adapt a concordance cosmology 
with (fl m ,n A ,h) = (0.3,0.7,0.7), and AB magnitudes 
with 2" diameter aperture, unless otherwise specified. 

2. OBSERVATION 

The new fully depleted red-sensitive CCDs were in- 
stalled in July 2008 on the Suprime-Cam. The quantum 
efficiency of the previous and new CCDs as well as trans- 
mission curve of the NB973 filter are shown in Figure 1. 
The quantum efficiency of the new CCDs at the central 
wavelength of NB973, 9755A, is ~ 0.6, about twice bet- 
ter than that of the previous CCDs. This substantially 
reduces the integration time required to reach a deeper 
limit than the SDF NB973 survey. 

We targeted another sky region studie d by the Subaru 
Obse rvatory Project, the SXDS field (jlurusaw a et al.1 
2008), because deep optical broadband BVRi'z' images 
taken by the Suprime-Cam are available. Also, deep nar- 
rowband NB816 (A c = 8160A and AA FW hm = 120A cor- 
responding to z ~ 5.7 Lya) and NB921 (A c = 9196A and 
AAfwhm =132A corresponding to z ~ 6.6 Lya ) filter 
image s of the SXDS were taken by lOuchi et al.l (2008, 
I2009al) . The entire SXDS field has an area of ~ 1.3 
square degree, consisting of five pointing of Suprime- 
Cam. They are named as SXDS-C, SXDS-N, SXDS- 
S, SXDS-E and SXDS-W with the central coordinates 
02 h 18 m 00. s 0, -05°00'00" (J2000) of the SXDS-C corre- 
sponding to the center of the entire SXDS field. More- 
over, this field has been observed with several differ- 
ent wavelength from X-ray to radio. Hence, multi- 
wavelength study of z = 7 LAEs is also possible. For 
example, the SXDS-C and most part of the SXDS-S 
have been obser ved with the UKIDSS U ltra Deep Survey 
(UKIDSS-UDS; ILawrence etaLl \2Wjb and the Spitzer 
legacy survey (SpUDS; PI: J. DunlopJEL However, the 
SXDS-C includes a few bright stars whose stellar ha- 
los contaminate large areas. Hence, to image the same 
area as the UKIDSS and the Spitzer surveys and to 
avoid those bright stars, we observed an area between the 
SXDS-C and SXDS-S, 02 h 18 m 00. s 0, -05°13'30" (J2000), 
as shown in Figure 2 by one pointing (34' x 27') of the 
Suprime-Cam. 

The NB973 imaging observations with the new red- 
sensitive CCDs were carried out at dark clear nights on 
25 and 26 October and 30 November in 2008. The sky 
conditions were photometric with the seeing of ~ 0."5- 
l."2. By taking the 30 minute exposure frames with 
8-point dithering patterns, we obtained a total of 17 
hours of imaging data. We have reduced the NB973 
imag e frames usi ng th e software SDFRED ([Ouchi et al.l 
l200l lYaei et al.1 [20 021 in the same standard manner as 
in Kashik awa et all (|2004l ) and lOta etaLl (|2008[ ). The 
dithered NB973 image frames were combined by dis- 
carding the frames with seeing size of > 1" because we 
confirmed that they degrade the final stacked image al- 
though they do not help improve the depth. The as- 
trometry of the combined image was then matched to 
those of the BVRi'z' SXDS-C and SXDS-S images to 
produce the two images we hereafter call NB973-SXDS- 
C and NB973-SXDS-S. The seeing size and integration 
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Fig. 2. — T he image of the Subar u/ XMM-Newton Deep Survey 
(SXDS) Field (Furusawa ct al. 2008). The north is up and the cast 
to the left. The SXDS consists of five pointing (totally ~ 1.3 square 
degree) of the Suprime-Cam named SXDS-C, SXDS-N, SXDS-S, 
SXDS-E and SXDS-W. We imaged the region between the SXDS-C 
and SXDS-S with the NB973 (the red squa re). The green, cyan and 
magenta lines show the UKlDSS-UDS (Lawrence ct al. 2003) and 
SpUDS fields (PI: J. Dunlop) with Spizer Space Telescope Infrared 
Array Camera 3.6^tm and 5.8/mi imaging and 4.5/xm and 8.0/im 
imaging, respectively. 

time of these final combined NB973 images are l."0 and 
13 hours. Spectrop hotometric standard stars LDS749B, 
GD248 and GD71 (jOkd 119901 ) were imaged during the 
observation to calibrate the photometric zeropoint of 
these stacked images, which was turned out to be NB973 
= 32.24 mag ADU -1 . The limiting magnitude reached 
NB973 < 25.4 at 5<7 with 13 hour integration, about 0.5 
magnitude deeper than our previous 15 hours of the SDF 
NB973 survey. 

3. PHOTOMETRIC CANDIDATES 

3.1. Photometry 

Using the NB973-SXDS-C and NB973-SXDS-S 
images, we made an NB973-detected object catalog 
and performed photometric analysis. Because the 
broadband BVRi'z' and narrowband NB816 and 
NB921 filter images of the SXDS-C and SXDS-S have 
smaller seeing sizes of ~ 0."7-0."8, we convolved 
these images to have the same seeing size of 1" as the 
NB973-SXDS-C and NB973-SXDS-S for the purpose 
of aperture photometry. The 2" aperture limiting 
magnitudes of these images at 3cr wit h the Galac- 
tic e xtinction corrected based on iSchlegel et al.1 
(p98h are (£, V, R, i', z', NB816, NB921) 
(28.00,27.71,27.52,27.58,26.54,26.52,26.37) for the 
SXDS-C and (28.24, 27.68, 27.62, 27.43, 26.36, 26.62, 26.37) 
for the SX DS-S. We used the SEx tractor software ver- 
sion 2.2.2 (B ertin fc Arnoutslll996h for source detection 
and photometry. The Suprime-Cam CCDs have a pixel 
size of ~ 0."2 pixel -1 . We considered an area larger 
than contiguous 5 pixels with a flux [mag arcsec ] 
greater than 2a to be an object. Object detection was 
first made in the NB973-SXDS-C and NB973-SXDS-S 
images and then photometry was done in the images of 
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Fig. 3.— The detection completenesses of t he NB973 SXDS (the 
present survey) and SDF (the previous survey IO"ta ct al. 20Q3) im- 
ages, calculated for every 0.5 mag bin. The completenesses do not 
reach 1.0 even for the objects with bright magnitudes. This is be- 
cause the blended or overlapped objects tend to be counted as one 
object by the SExtractor. The completenesses are corrected when 
the number and luminosity densities of z = 7 LAEs in the SXDS 
and SDF are calculated in §4. Our current survey in the SXDS is 
much deeper and more complete than our previous survey in the 
SDF. 

the other wavebands for SXDS-C and SXDS-S, respec- 
tively, using the double-image mode. We measured 2" 
diameter aperture magnitudes of detected objects with 
MAG-APER parameter and total magnitudes with 
MAG_AUTO. The final object catalog was constructed 
by combining the photometry from both SXDS-C and 
SXDS-S, detecting the total of 230,044 objects down to 
NB973 < 25.4 {5a). 

3.2. Selection of Candidate z — 7 LAEs 

We used the same color criteria as ones adopted for our 
previous NB973 survey in the SDF (jOta et al.ll2008h to 
select out candidate z = 7 LAEs in the SXDS. Namely, 

(1) B,F,.R,i / ,NB816,z / ,NB921 < 3ct; NB973 < 25.4 

(2) B, V, R, NB816, NB921 < 3cr; %' - z' > 1.3; 
z' - NB973 > 1.0; NB973 < 25.4 

We found 1274 and 18 objects satisfying the criterion (1) 
and (2), respectively. 

However, most of them are located in the noisy rea- 
gions near the edges of the NB973-SXDS-C and NB973- 
SXDS-S images where the signal-to-noise (S/N) ratio is 
low. This implies that most of them could be noises that 
appear only in the NB973 images and not in the other 
wavebands. To examine if they are noises, we created the 
negative NB973-SXDS-C and NB973-SXDS-S images by 
multiplying each pixel value by —1 and performed source 
detection using SExtractor with the detection and analy- 
sis thresholds of 5a, instead of 2a, because our candidate 
selection adopted NB973 < 25.4 (5a) and the sky back- 
ground rms's do not change even if we change the positive 
images into negative ones. The top and right edges of the 
NB973-SXDS-C image and the bottom and right edges 
of the NB973-SXDS-S images were dominated with nega- 
tive > 5a detections, which are considered noises. These 
edge regions coincide with the locations where most of 
the sources selected with the z = 7 LAE candidate cri- 
teria (1) and (2) above distribute. Hence, we trimmed 
these edge regions off the NB973-SXDS-C and -S im- 
ages. After that, we also measured the limiting magni- 
tudes near the edge of the images and found that the 
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Fig. 4. — z' — NB973 (2" aperture mags.) color as a function 
of NB973 (total) magnitude of the objects detected in our survey 
field (shown by dots). The dotted curve shows 3a error track of 
z'— NB973 color. The horizontal solid line is a part of our color 
selection criterion, z'— NB973 > 1.0. The vertical dashed line 
indicates the detection limiting magnitude of our survey, NB973 
= 25.4(5cr). The diagonal dashed line is the 2a limits of z'-band 
2" aperture magnitude. The three z = 7 LAE candidates are 
denoted by the filled circles with the upper arrows showing the 2a 
limits on z'— NB973 colors. 

regions within 1000 pixels (~ 3/4) from right edges of 
the NB973-SXDS-C and -S images are 0.3 mag shallower. 
Thus, we also removed these region from the images. 

Consequently, we were left with three objects located 
in NB973-SXDS-S image, all of which satisfy the can- 
didate selection criterion (1). We visually inspected all 
the broadband and narrowband images of each of these 
three objects to remove obviously spurious sources such 
as columns of bad pixels, pixels saturated with bright 
stars, noise events of deformed shapes and scattering pix- 
els having anomalously large fluxes. None of the these 
applied to the three objects. We also examined if the ob- 
jects are not seen in the BVRi', NB816 and NB921 be- 
cause such objects are considered low-z interlopers since 
flux at the wavelength shortward z = 7 Lya should be 
mostly absorbed by the IGM. None of the three were seen 
in any of the wavebands blueward z — 7 Lya. 

However, because these selected objects are only de- 
tected in the NB973 image and not in the other wave- 
bands, they might include transient or variable objects 
such as active galactic nuclei (AGNs) and supernovae 
that were too faint to be detected at the time of BVRi' z' , 
NB816 and NB921 imaging observations of the SXDS 
but became brighter than 5a at the time of our NB973 
imaging. Also, the three objects might still include noise 
events that exist only in the NB973 image and that re- 
semble real objects in shape. The color selection and 
the visual inspection above alone cannot remove these 
contaminants. 

To reduce the number of these contaminants as much 
as possible, we conducted another type of photometric 
analysis as follows. We created three additional NB973 
SXDS images: (1) a stacking of NB973 image frames 
(11 hours in total and PSF of l."2) taken on 25 and 26 
October 2008, (2) a stacking of NB973 image frames (6 
hours in total and PSF of l."2) taken on 30 November 
2008 and (3) a stacking of all the NB973 image frames 
including the ones with PSF > 1" (17 hours in total and 
PSF of l."2) taken in both October and November 2008. 
Then, we measured the 2" aperture fluxes of the three 
objects in each of the images (l)-(3). We found that each 
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Fig. 5. — The multi-waveband images of the three candidate z = 7 LAEs. The size of each image is 20" X 20". They are detected only 
in NB973. 



object was detected in all the images (l)-(3) at 4. 3-5. Oct, 
5.0-7.8cr and 3.0-4.4<r significances, respectively. On the 
other hand, we also created the negative images from 
the images (f)-(3), performed source detection in the 
same way above and found that negative > 5<r noises 
were not detected in all the images (l)-(3) in the same 
locations. Hence, the three objects are not likely to be 
noises. Meanwhile, the difference in NB973 flux between 
the images (1) and (2) (i.e., between the October and 
November 2008) in each of the three objects was 0.32er, 
1.7<7 and 0.66er, respectively. Since the flux differences 
are within ±lcr photometric errors, all the three objects 
did not show any significant variability in one month. 
Hence, they might not be variable objects. However, we 
cannot completely exclude the possibility of them being 
variable objects because their time variability might be 
longer than one month. 

Eventually, we were left with three candidates and 
this is an upper limit on the number of z — 7 LAEs 
in the NB973 SXDS field since they might still include 
variable objects with slow variability. In ou r previous 
NB973 survey in the SDF, lOta et al.1 (|2008f ) estimated 
that the number of variables could be ~ 3.8-4.5 in the 
876 arcmin 2 of sky area. If this number is correct and 
also applies to the NB973-SXDS-C and NB973-SXDS- 
S images (totally 824 arcmin 2 , see below), all the three 
objects selected in the SXDS could be variable objects. 
However, to securely identify them, the follow-up spec- 
troscopy is of course required. The color-magnitude di- 
agram (V- NB973 vs. NB973) of the three LAE candi- 
dates and all the objects detected down to NB973 = 25.4 
is plotted in Figure 4. The images of the three LAE can- 
didates and their photometric properties are shown in 
Figure 5 and Table [TJ respectively. 

To remove the noises in the process of the candidate 
selection, we trimmed the low S/N edge regions off the 
NB973-SXDS-C and NB973-SXDS-S images. After this, 
the total area of the NB973-SXDS-C and NB973-SXDS- 
S images became 824 arcmin 2 . The comoving distance 
along the line-of-sight corresponding to the redshift range 
6.94 < z < 7.11 for LAEs covered by NB973 filter is 58 
Mpc. Therefore, we have probed a total of 3.0 x 10 5 
Mpc 3 volume in our survey. 

3.3. The Detection Completeness 

The number of objects we can detect decreases as the 
fluxes of objects become fainter. To examine what frac- 



tion of objects in the NB973 image the SExtractor can 
detect or fails to detect down to the limiting magnitude of 
NB973 < 25.4, we measured the detection completeness 
of our photometry for the NB973-SX DS-C and NB973 - 
SXDS-S images in the same way as IQta et all ((2008). 
The detection completeness is the ratio of the number of 
detected objects to that of objects actually imaged in the 
NB973-SXDS-C and NB973-SXDS-S images. The result 
is shown in Figure |3l The completeness is ~ 72% at our 
detection limit of NB973 = 25.4. This completeness was 
corrected when the number and luminosity densities and 
Lya LF of z = 7 LAEs were calculated in the §4. Fig- 
ure [3] also shows the completeness of our previous NB973 
survey in the SDF for comparison. It is clear that the 
present survey in the SXDS is much deeper and more 
complete. 

3.4. Lya Luminosity and Star Formation Rate 

To derive the Lya LF of z = 7 LAEs, we need to know 
their Lya fluxes. Now that we know the NB973 magni- 
tudes of our z = 7 LAE candidates, we can estimate their 
Lya fluxes, -F(Lya), Lya luminosities, L(Lya), and cor- 
responding star formation rates (SFRs), SFR(Lya), by 
directly converting the NB973 magnitudes into F(Lya). 
Note that these quantities could be upper limits because 
the narrowband NB973 filter might possibly detect not 
only the Lya emission of a LAE but also its UV contin- 
uum flux if it is strong. In our present survery, we have 
neither sufficiently deep near-infrared images nor spec- 
troscopic data for our three z = 7 LAEs for the estima- 
tion of UV continuum fluxes to correct the narrowband 
fluxes for the continuum contribution. The -F(Lya), 
L(Lya) and SFR(Lya) calculated here are listed in Table 
El We convert F(Lya) into L(Lya) by assuming that the 
redshift of Lya emission is z = 7.02 corresponding to the 
central wavelength of NB973 (9755A). The uncertainty 
of L(Lya) due to the location of Lya emission (z — 6.94- 
7.10) within the NB973 passband (A = 9655-9855 A) is 
AlogL(Lya) < 0.02. To estimate SFR(Lya), we used 
the following rela tion derived from Kennicutt's equation 
( Kennicutt]ll998f) w ith the case B recombination theory 
(|Brocklehurstlll971h . 

SFR(Lya) = 9.1 x 10" 43 i(Lya)M Q yr" 1 (1) 

It should be noted that this relation holds good in the 
case where Lya emission is not attenuated by neutral 
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TABLE 1 

Photometric properties of candidate z = 7 Lyq emitters 



Object 


R.A. (J2000.0) 


Decl. (J2000.0) 


i> 


NB816 


z> 


NB921 


NB973 


NB973 (total) 


NB973-SXDS-S-5729 


02:17:48.97 


-05:27:18.56 


> 27.87 


> 27.06 


> 26.80 


> 26.81 


25.24 


24.52 


NB973-SXDS-S-113268 


02:17:59.54 


-05:14:07.61 


> 27.87 


> 27.06 


> 26.80 


> 26.81 


24.84 


24.57 


NB973-SXDS-S-66924 


02:17:57.85 


-05:18:47.46 


> 27.87 


> 27.06 


> 26.80 


> 26.81 


25.30 


25.10 



NOTE: Units of coordinate are hours: minutes: seconds (right ascension) and degrees: arcminutes: arcscconds (declination) using J2000.0 
equinox, i' , NB816, z' , NB921 and NB973 are all 2" aperture magnitudes while NB973 (total) is a total magnitude. Magnitudes are 
replaced by their 2a limits if they are fainter than the limits. 



IGM. Lya emission could be attenuated by neutral IGM 
at z > 6. Hence, we basically discuss L(Lya) and its 
volume density pLya as physical properties of z = 7 LAEs 
in the subsequent sections and list SFR(Lya) and its 
volume density SFRDL yQ just for reference. 

4. DISCUSSION 

In our previous survey, we detected one z ~ 7 LAE 
to the limit of L(Lya) > 1 x 10 43 erg s" 1 in the SDF. 
The number and Lya luminosity densities of z = 7 LAE 
were significantly smaller than those of LAEs at z = 5.7 
and z — 6.6, suggesting that neutral fraction might be 
larger at higher redshifts beyond z ~ 6. However, the 
SDF survey had three weaknesses: shallow de pth, small 
statis tics and observation of only one sky field f Ota et al.l 
2008). We have now detected three LAE candidates in a 
different sky field, SXDS, to 0.5 mag deeper limit than 
the SDF survey. Hence, we derive a deeper Lya LF of 
z = 7 LAEs with improved statistics to investigate what 
it implies for galaxy evolution and reionization and to 
sec if it supports our previous result or how the result 
changes in the different sky field. 

4.1. Lya Luminosity Function of z — 7 LAEs 

In H'SAl we estimated the upper limit on L(Lya) for 
each LAE candidate. Based on this, we now derive the 
upper limit on Lya LF of z = 7 LAEs at z — 7 to our 
survey limit of NB973 < 25.4 corresponding to L(Lya) > 
9.2 x 10 42 erg s _1 (also assuming that Lya emission is 
at z = 7.02 corresponding to the central wavelength of 
NB973 filter 9755 A). We estimate the number density of 
LAEs by simply dividing the cumulative number of LAEs 
in each L(Lya) bin by the survey volume of 3.0 x 10 5 
Mpc 3 estimated in §2. For the error of number density in 
each bin, we include Poisson err ors and co s mic y ariance 
estimated in the same way as in IQta et all d200p) . That 
is, we used the Tables 1 and 2 in IGehrelsl (|1986l ) for the 
estimate of Poisson erro rs and the semi-analytic model of 
iSomerville et al.l (|2004D for the cosmic variance. Also, we 
correct the number density and its error for the detection 
completeness estimated in §2 and Figure 3. The derived 
upper limit Lya LF of z = 7 LAEs in the SXDS as well 
as LFs of LAEs at z = 5.7 6.6 and 7.0 in the SDF o r 
SXDS fromlShimasaku efah| (l2006h.lOuchi et al.l (12001. 
Kashikawa ct al 



(2006) and 



Ota et al.l (|2008f) are shown 



in Figure [5] It should be noted that while our LF is the 
upper limit without correction for contribution from UV 
continuum flux, all the other LFs are corrected for this 
using either the deep broadband images redward Lya or 
spectroscopy data (See each reference). 

The upper limit z = 7 LF is similar to z = 5.7 LFs 
at its brighter end but shows deficit from z = 5.7 at the 
fainter luminosity. The fainter part of the upper limit z = 
7 LF is similar to the z — 6.6 LF. The number density 



of the SXDS upper limit z — 7 LF is larger than that of 
the SDF z — 7 LF based on the spectroscopy. However, 
Figure [6] also shows that if we assume that .F(NB973) = 
i^(Lya) for the SDF z = 7 LF, it is similar to the SXDS 
z = 7 LF. This suggests that the field-to-field variance 
in z = 7 Lya LF might not be so significant between 
the SDF and SXDS. However, to investigate the field-to- 
field variance more precisely, we should know accurate 
Lya luminosities of the three z = 7 LAE candidates with 
spectroscopy. 

Figure [7J and Table [3] show the LAE number den- 
sity riL yQ , and Lya luminosity density pLya estimated 
from our three z = 7 LAE candidates to the survey 
limit of L(Lya) = 9.2 x 10 42 erg s" 1 . They are 
nlf a ~ 1.9±g;g x 10" 5 Mpc" 3 and plf a ~ 1.9±g;| x 10 38 
erg s" 1 Mpc" 3 , respectively. Also, in Figure [7] and Table 
[3] to the same limit, the number and Lya luminosity 
densities at z = 5.7 and z — 6.6 are estimated to 

\)q ( <n z ~ 5.7phot ^2— 5.7phot ^z— 6.6phot ^z— 6.6spec 



/ 2— 5.7phot 
V n Lya ' 
d+4.2 fi o + 1.5 2 r + 1.6 n i+i.M 
5 -3.6' °-°-1.5> z '°-1.4' 1 - L -0.8J X 

/ z — 5.7phot z— 5.7phot z— 6.6phot z— 6.6spcc\ 
\Phya ' Phya ! Phya > ^Lya I 

5.9 n n+1-9 o 0+1-9 i K+l^ ^ t n 38 „„„ e ,-l 



. 5.7phot z— 6.6phot 
n Lya ' ""Lya ' n Lya 

+ 1.^ . 10 -5 



Mpc 



-3 



and 

(12.4ti;;,9.0ti;g, 2.9113, 1.5li;§) x 10 38 erg s" 1 Mpc" 
from the z = 5.7 Lya LFs based on photometric LAE 
sampl es in SDF and SXDS de rived by iShimasaku et al.l 
(2006) and lOuchi et"aH ([2001 and the z = 6.6 Lya LFs 
based on photometric and s pectro scopic LAE samples 
derived by Kashik awa et al.l ()2006h . respectively. These 
densities were ca lculated by integrating the best-fit 
iSchechteFI (|1976f ) LFs derived in these studies and 
each uncertainty consists of Poisson error an d cosmic 
varian ce. Again, we used the Tables 1 and 2 in iGehrels 1 
(1986) for the estima te of Poisson e r rors and the 
semi-analytic model of ISomerville et al.l (|2004l ) for the 
cosmic variance. Similarly, densities at 3 < z < 5.7 were 
also calculated using the best-fit Schechter Lya LFs 
from previous studies (|Ouchi et al.l 120081 : iDawson et al.l 
120071) and shown in Figure [Jj The LAE number 
and Lya luminosity densities do not change much at 
z < 5.7 but do decrease at z > 5.7. The upper limit 
on number and Lya luminosity densities at z = 7 are 
ntfjnt; 5 * 7 - 7.7-54% and pifjpl= 5 a 7 ~ 5.5-39% of 
those at z = 5.7. 

We have discussed differences in n^ ya and ph ya be- 
tween z = 7 and later epochs based on the upper limit 
z = 7 Lya LF in the SXDS. This LF was derived from 
Lya fluxes of z — 7 LAE candidates directly converted 
from NB973 magnitudes without correcting for UV con- 
tinuum fluxes. Even so, the z = 7 LF shows deficit 
in number and Lya luminosity densities from z — 5.7. 
Hence, in reality, if the UV continuum fluxes are cor- 
rected, the z = 7 LF might show more deficit. 

4.2. Comparison with a Galaxy Evolution Model 
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TABLE 2 

Luminosities and SFRs of candidate z = 7 Ly« emitters 

Object F(Lya) a L(Lya) b SFR(Lya) c 

(10~ 17 erg s" 1 cm" 2 ) (10 42 erg s" 1 ) (A^yr" 1 ) 

NB973-SXDS-S-5729 3^6 20.7 ±0.6 18.8 ±0.5 

NB973-SXDS-S-113268 3.4 19.7 ±0.5 18.0 ±0.5 

NB973-SXDS-S-66924 21 12.1 ±0.3 11.0 ±0.3 

"The Lya flux calculated assuming that F(Lya) = F(NB973). 

'The Lya luminosity converted from F(Lya) assuming that Lya emission is located at the center of NB973 passband (i.e., z = 7.02). The 
uncertainties correspond to the other locations of Lya emission (z = 6.94-7.10) within the NB973 passband (A = 9655-9855A). 
c The SFR corresponding to L(Lya). 



The deficit in densities from z = 5.7 to z = 7 
could reflect both or either of galaxy evolution between 
these epochs and attenuation of Lya emission by neu- 
tral hydrogen remaining at the epoch of reionization. 
It is possible to estimate how much the number den- 
sity or Lya luminosity density changes due to galaxy 
evolution by using empirically-calibrated semi-analytic 
galaxy evolution models. In Figure HJ we also plot 
Lya LFs of LAEs at z = 5.7, 6.6 and 7 expected 
in the case of the IGM transmission for Lya photons 
^Lya 1 = 1 ( or equivalently neutral fraction xhi = 0), 
w hich we calculated by using a LAE evo l ution model 
of iKobavashi. Totani fc Nagashimal (120071 120101 here- 
after KTN10 model) . iKobavashi. Totani fc Nagashima I 
(|2007l I2010D constructed this model by incorporating 
new modeling for an escape fraction of Lya photons from 
galaxies into a recent hier archical clustering model of 
iNagashima fc Y oshii (2004), physically considering dust 
extinction of Lya photons and the effect of galaxy-scale 
outflows. The KTN10 model was empirically calibrated 
to fit the observed Lya LF of z = 5.7 LAEs in the SDF 
(jShimasaku et all [2006). It is worth noting that, with 
the consistent set of model parameters, the model nat- 
urally reproduces all of the available observed data of 
the LAEs (i.e., Lya LF, UV LF, and equivalent width 
distribution) in the redshift range of z ~ 3-6 under the 
standard scenario of hierarchical galaxy formation. 

The KTN10 model z ~ 7 LF shows lower number 
density than z = 5.7 LFs due to the galaxy evolution. 
The observed upper limit z — 7 LF is consistent with 
the KTN10 z = 7 LF at the bright end but shows a 
deficit at the fainter end. Integrating the KTN10 LF 
to our survey limit of LiLya) > 9.2 x 10 42 erg s , 
we obtain the number and Lya luminosity densities of 
<;I' KTN1 ° - 7.2 x 10^ Mpc-3 and p^'™ ~ 
9.5 x 10 38 erg s _1 Mpc -3 as also shown in Figure [7] and 
Table [3J Thus, the observed number and Lya luminos- 
ity densities are factors of nf^l/n^~^' KTN10 ~ 0.14-0.39 

and Piya/Ph^a KTN1 ° ~ 0.17-0.36 lower than the model- 
predicted densities of z — 7 LAEs in the environment 
where shi = 0. 

This density discrepancies between the model and ob- 
served LFs might possibly reflect the attenuation of Lya 
luminosity by neutral hydrogen. If we assume that this 
is the case, we can estimate the IGM transmission for 
Lya photons T L ™ at z = 7 with T L ™ = observed 
L(Lya)/KTN10 model L(Lya). We shifted the KTN10 



z = 7 LF by attenuating its Lya luminosity L(Lya) and 
compared it with our observed z = 7 LF to find the 
possible range of T^^ 1 . In Figure [HI we also show the 
attenuated KTN10 z = 7 LF, which still passes through 
the observed z — 7 LF data within error at the faintest 
end and gives the maximum possible attenuation factors 
of T L y^ = 0.4. Hence, the Lya luminosity is attenuated 

by a factor of 0.4 < T[™ < 1 (see also Table©. 

4.3. Implication for Reionization 

In £14.21 we estimated the effect of galaxy evolution on 
the decrease in the number density by using the KTN10 
model. Even after correctiong for the galaxy evolution 
effect, the number density of z = 7 LAE Lya LF still 
shows a deficit. This could be due to Lya photon absorp- 
tion by neutral IGM, implying a possible evolution of the 
IGM neutral fraction between z = 5.7 and 2 = 7. If we 
can translate T^^f obtained in £ 14.21 to neutral fraction 
at z = 7, we might be able to obtain some quantitative 
implications for reionization. However, this conversion 
process is highly complicated and uncertainty is large be- 
cause we have to rely on currently proposed models that 
quant ify Lya emis s ion line attenuation by neutral IGM 
(e.g lSantoi 120041: IDiikstra et all 120071 IMcQuinn et~aTl 
2007). As we did in IQta et all (|2008[ ) and for compari- 
son with their result from the z = 7 LAE survey in the 
SDF, we use the dynamical model with a reasonable ve- 
locity shift of the Lya line by 3 60 km s" 1 redward of the 
systemic velocity (|Santosll2004D . which provides as 
a function of neutral fraction x-gj.- 

In $L2l we obtained T[™ = 0.4-1 from the com- 
parison of the KTN10 model LF and the upper limit 
on z — 7 Lya LF in Figure [51 If we appl y this 
T L y^ to the Lya flux attenuation model of ISantosI 
(2004), the neutral fraction at z = 7 would be < 
x m 7 S 0.63. This range of x^ 7 is consistent with 
xgf 7 ~ 0.32-0.62 obtained by ou r previous NB973 sur- 
vey in the SDF (IQta et all 120081) whic h also used the 
IKobavashi. Totani fc Nagashimal (|2007t) LAE evolution 
model, the previous version of the KTN10 model. Since 
the difference i n predicted number density of z = 7 Lya 
LF between the lKobavashi. Totani fc Nagashi ma I (|2007l ) 
model and the KTN10 model is negligibly small to our 
survey limit, the consistency of the result between the 
SXDS and SDF implies that the field-to-field variance 
in reionization state is not so significant between these 
different sky fields of comparable areas. 

Meanwhile, the Gunn-Peterson troughs 
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z=5.7 model 

z=6.6 model 

z=7.0 model 

z=7.0 model 0.40 x L(Lya) 

Z=5.7 SXDS PtlOt : :! ' 

z=57SDFPhot ^-U- J . 

z=6.6 SDF Phot : • : 
z=6.6 SDF Spec a 
z=7 SDF Spec 
This study SXDS z=7 Phot f 
z-7 SDF Phot Q 




Age of the Universe (Gyr) 
2.03 1.66 1.31 0.97 



42 42.5 43 43.5 

log[L(Lya) (erg s" 1 )] 

Fig. 6.— Cumulative Lya LFs of z = 5.7 LAEs in the SDF 
UShima saku ct al. 2006]) and S XDS IjOuchi et al.l I200SI '). z = 6.6 
LAEs Eyhikawa et alJ EEl) and z = 7 LAEs in the SDF 
llOta et aU2008l ) and SXDS (this study based on the 3 candidates). 
The LFs based on the photmetry and spectroscopy samples are la- 
beled as Phot and Spec. The z = 7 SXDS LF is the upper limit LF 
based on the assumption that F(hya) = .F(NB973) in each 2 = 7 
LAE candidate while L(Lya)'s of all the other LFs are corrected for 
the contribution from UV continuum fluxes (See esch reference). 
There are 4 data points for our 3 z = 7 candidates because the last 
bin corresponds to our survey limit of Li; m ; t (Lya) = 9.2 X 10 42 erg 
s -1 (shown by the arrow) and the faintest candidate is brighter 
than this li mit. A ll the errors include cosmic variance and Poisson 
errors (See N4.1H . All the data and errors are corrected for their 
detection completeness (See Figure |3). The long-dashed, short- 
dashed and solid curves are the intrinsic (i.e., T^^f = 1) Lya LFs 
at z = 5.7, 6.6 and 7, resp ectively, predicted by the KTN10 LAE 
evolution model (See i|4.2H . The thin solid curve is the KTN10 
model 2 = 7 LF with its Lya luminosity attenuated by a factor of 
0.4 to be consistent with the data points with its errors that result 
in the largest attenuation of the LF. 

(|Gunn fc Peterson! I1965D found in z ~ 6 quasar 
spectra suggest that reionization ended at z ~ 6 with an 
estimated neutral fraction at x 6.2, ^j^i 

6 2 - 0.01-0.04 

(|Fan et al.1 120061 ). Also, spectral analyses of z ~ 6.3 
and 6.7 GRBs imply that reionization is not yet com- 
plete at z ~ 6.3 and 6.7 with x^ e 3 < 0.17-0.6 and 
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xgf- ' > 0.35 (jTotani et all [2001 iGreiner et alJl2009h . 
These results, together with the z = 7 neutral fraction 
we derived, imply that xm could be higher at the 
earlier epoch beyond z ~ 6 at which reionization is 
believed to have ended. However, it should be again 
noted that our quantitative results from LAE surveys 
are model-dependent and thus considered one possible 
interpretation. All the densities, IGM transmission and 
neutral fraction obtained in i j4.lH4.3l (nL yQ , Phya , 
=7S| are summarized in Table [3J 



T-rIGM 
Lya < 



5. CONCLUSION 



We conducted a deep narrowband NB973 survey of 
z = 7 LAEs in the SXDS using the fully depleted CCDs 
newly installed on the Subaru Telescope Suprime-Cam, 
which doubled the detector sensitivity to z = 7 Lya emis- 
sion at ~ 1/im. Integrating 13 hours, we reached a survey 
limit of NB973 = 25.4 (5a), 0.5 mag deeper than our pre- 
vious survey in the SDF. To this limit, we detected three 
candidate z = 7 LAEs, all of which are detected in all 
the multi-epoch NB973 imaging data and thus probably 
real objects. 

With the survey depth and the sample statistics sig- 
nificantly improved compared with our previous NB973 
survey in the SDF, we derived deeper but upper limit 
Lya LF at z = 7, assuming that F(NB973) = F(Lya). 
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Fig. 7. — The number density WLyai Lya luminosity density 
Phya an d star formation rate density SFRDL yQ of LAEs at 2 = 
5.7, 6.6 and 7 derived from the current study and the latest 
Subaru/Suprime-Cam LAE surveys in the SDF and SXDS and 
those at 3 < 2 < 5.7 from literature to ^limit 

(Lya) = 9.2 X 10 42 

erg s _1 . The densities at 2 = 7 (large filled circles) are calcu- 
lated from our 3 2 = 7 LAE candidates in the SXDS with the 
assumption that F(Lya) = F(NB973). The densities at all the 
other redshifts are calculated by integrating the the best-fit Lya 
Schcchter LFs derived f rom the 2 = 6.6 SDF photometric and 
spectroscopic samples of Kashikawa ct al. (2006, open circles and 
triangles, respect i vely), the 2 = 5.7 SDF photometric sample of 
IShimasa ku ct al. (200 6], open s quare s), the 2 = 3.1, 3.7 and 5.7 
SXDS sa mples of Ouchi ct al. ( 2008,, diamonds) and the 2 ~ 4.5 
sample of Dawson et al. (2007, pentagon). All the parameters of 
the Schechter LFs (normalization <j>* , characteristic luminosity L* 
and faint end slope a) used for the calculations of the densities 
shown here are the same as those used for the similar density cal- 
culations in Ota ct al. (2008). Each horizontal error bar shows the 
redshift range of each survey. The vertica l err or bars include both 
cosmic variance and Poisson errors (See Sj4Tj. The plus symbols 
at 2 = 5.7, 6.6 and 7 with solid lines show the expected densi- 
ties obtained by integrating the intrinsic Lya LFs predicted by the 
KTN10 model. At 2 > 5.7, the densities decrease and smaller than 
the model-predicted values, implying that the Lya lines might be 
attenuated by the possibly increasing neutral IGM at reionization 
epoch. 

To our survey limit L(Lya) = 9.2 x 10 42 erg s~ 1 (corre- 
sponding to NB973 = 25.4), the z = 7 LF is similar to 
the z — 5.7 LAE Lya LF at the brighter end but shows 
a clear deficit from the z = 5.7 at the fainter end. The 
fainter end of the z = 7 LF also looks similar to the 
z — 6.6 LAE Lya LF. The number and Lya luminos- 
ity densities at z = 7 are nf= 7 ~ l^t° i X 10~ 5 Mpc~ 3 



and pl y l 



i.9±8:l x 10 38 



Lya 

erg s 



-0.9 
-1 A/T^„-3 



Mpc" 



ry-iZ 7 / ry-. Z 5.7 

«t ya/'tLya 



Lya/ "Lya ~ 7.7"54% and P^J 



They are 
5.5-39% of 

the densities at z = 5.7. 

The decrease in number and Lya luminosity densities 
from z = 5.7toz = 7 could be due to either or both of 
galaxy evolution and suppression of Lya line flux by neu- 
tral hydrogen. Our observed z — 7 Lya LF shows deficit 
from the predicted z — 7 LF corrected for the possible 
density decrease due to galaxy evolution calculated by 
using a LAE evolution model (the KTN10 model). If 
we attribute the deficit to the attenuation of Lya flux by 
neutral hydrogen, the IGM transmission for Lya photons 
at z = 7 would be 0.4 < T[™ < 1. 

Our current result as to strengthens the possibil- 

ity that reionization might not be complete yet at z — 7 
and supports the previous results from the z ~ 6 quasars 
and z > 6 GRBs, which suggest that reionization might 
have completed at z ~ 6. The Lya LF, trend of decrease 
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TABLE 3 

LAE Number and Lyo Luminosity Densities, IGM Transmission and Neutral Fraction 

Lya Luminosity Function" ^ ^ T™~ c xgf d 

(10~ 5 Mpc~ 3 ) (10 38 erg s" 1 Mpc~ 3 ) 

^ - 7.0 Upper 1.9± u % 1.9+^ = = 



5.7 SDF Phot 8.8±|'| 12A -ii 

+i's an+i-9 



z = 5.7 SXDS Phot 6.8IJ;J 9-0-i _ 

^ = 6.6 SDF Phot 2.5±i;| 2.9±^;§ 

z = 6.6 SDF Spec 1-1±08 

^ = 7.0KTN10 7.2' 9.5' 0.40-1.00 0.00-0.63 



"These LFs correspond to the ones in Figure [6] 

''Number and Lya luminosity densities calculated to our survey limit of L(Lya) = 9.2 X 10 42 erg s — 1 (See £|4,ll and l4.2ll , 
C IGM transmission for Lya photons at z = 7 derived from the comparison of z = 7 SXDS LF and 2 = 7 KTN10 LF (Sec 34.21 1. 
d Neutral fraction at z = 7 obtained by appying the T^M to lSantosl 1200411 model. 



ill TlLya and /^Lya 

from z = 5.7, and T[f™ at z = 7 
investigated in the present deep NB973 survey in the 
SXDS are all consistent with those studied in our previ- 
ous NB973 survey in the SDF. This implies that field-to- 
field variance in these properties as well as reionization 
status between the two different sky fields of comparable 
areas, the SXDS and SDF (824 and 876 arcmin 2 , respec- 
tively, corresponding to the survey volume of ~ 3 x 10 
Mpc 3 ), is not significant. However, observing two sky 
fields is statistically small and increasing the number of 
survey fields is still important to investigating the field 
variance of reionization state more accurately. 

The results regarding z — 7 LAEs presented here are 
based on deeper survey limit and larger sample statistics 
than the SDF and our current conclusion is robuster. 
However, the sample is constructed from the photomet- 
ric candidates, and the follow-up spectroscopy of them 
are indispensable to know how many of them are real, to 
estimate Lya flux and luminosity precisely and thus to 
increase the accuracy of the study on the galaxy evolu- 
tion and reionization. Even so, our current study shows 
that even if all the three candidates would be spectro- 
scopically identified to be real LAEs at z = 7, the fact 
that the number and Lya luminosity densities decrease 
from z ~ 6 to z = 7 will not change. If some of the can- 
didates turn out to be unreal, the densities will be even 
smaller. Hence, something could be really going on be- 



yond z ~ 6, which might possibly be reionization. This 
will still remain an open question for the time being, be- 
cause our constraint on reionization at z = 7 depends 
on one or two models: a LAE evolution model and/or 
a Lya flux attenuation model. The improvement of the 
accuracy of these semi-analytical and theoretical models 
in the future will help answer the question. 
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